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Abstract. The IRAM Interferometer has been used 
to detect the submillimeter source HDF 850.1 found 
by Hughes et al. (1998) in the Hubble Deep Field. 
The flux density measured at 1.3 mm (236 GHz) is 
2.2±0.3(ler) mjy, in agreement with the flux density mea- 
sured at the JCMT. The flux densities and upper limits 
measured at 3.4, 2.8, 1.3, 0.85, and 0.45 mm show that the 
emission comes from dust. We suggest that the 1.3 mm 
dust source is associated with the optical arc-like feature, 
3-593.0, that has a photometric redshift z w 1.7. If the 
HDF 850.1 is at this redshift and unlensed, its spectral en- 
ergy distribution, combined with that of 3-593.0, matches 
closely that of the ultraluminous galaxy VII Zw 31. An- 
other possibility is that the dust source may be gravita- 
tionally lensed by the elliptical galaxy 3-586.0 at z w 1. 

The position of the dust source agrees within the errors 
with that of the tentative VLA radio source 3651+1226. 

Key words: galaxies: structure - galaxies: individual 
(HDF-850.1) - cosmology - galaxies: ISM 
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1. Introduction 

Hughes et al. (1998) observed the Hubble Deep Field with 
the SCUBA multi-beam bolometer array at the James 
Clerk Maxwell Telescope (JCMT) at 850 and 450 ^m. 
Apart from HDF 850.4 (and possibly HDF 850.3) the 
SCUBA survey did not detect the mid-IR sources found 
in the Hubble Deep Field with the ISO satellite (Rowan- 
Robinson et al. 1997; Mann et al. 1997; Aussel et al. 1999). 
In general these ISOCAM detections are associated with 
z < 1 galaxies, while the submillimeter survey appears 
to have revealed a higher-redshift population of dust en- 
shrouded starburst galaxies. Further SCUBA surveys of 
other fields (e.g., Barger et al. 1999a, b; Eales et al. 1999; 
Lilly et al. 1999) indicate however that many of the sub- 
mm sources have redshifts smaller than 1, and that much 
of the star-forming activity in galaxies has occurred rela- 
tively recently, at z ~ 2. 

Richards (1999) argued against high redshifts (z > 2) 
for the HDF sub-mm sources, from plausible associations 
of the SCUBA sources with nearby radio sources detected 
in his VLA survey at 20 cm (1.4 GHz) and the identifica- 
tion of some of these non-thermal radio continuum sources 
with optical galaxies in the Hubble Deep Field. From the 
radio luminosity function of starburst galaxies, Richards 
argued that if the submm dust sources had cm-radio detec- 



2 



D. Downes et al.: Proposed identification of Hubble Deep Field submillimeter source HDF 850.1 



Table 1. Source Positions 



R.A. 


Dec. 




A6> from 


Redshift 


Remark 


Refs. 


Source (J2000) 


(J2000) 


HDF850.1 


2 






Proposed identification: 














HDF-850.1 12 h 36 m 51.98 s ± 0.04 s 


62°12'25.7" 


± 0.3" 


— 





1.3 mm dust source 


1 


3-593.0 12 h 36 m 51.86 s ± 0.06 s 


62°12'25.6" 


± 0.4" 


0.8" 


1.73 to 1.76 


"arc 1" 


2, 6, 7 


3-593.2 12 h 36 m 51.94 s ± 0.06 s 


62°12'25.0" 


± 0.4" 


0.8" 


1.82 


"arc 2" 


2, 6 


3651+1226 12 h 36 m 51.96 s ± 0.14 s 


62°12'26.1" 


± 1.0" 


0.4" 


— 


VLA (3 cm) ? 


3 


Other nearby optical sources within 3": 












3-586.0 12 h 36 m 52.10 s ± 0.06 s 


62°12'26.3" 


± 0.4" 


1.0" 


1.0 to 1.2 


elliptical 


2, 6, 7 


3-577.0 12 h 36 m 52.26 s ± 0.06 s 


62°12'27.2" 


± 0.4" 


2.5" 


2.88, 2.89, 3.36 ? 


"counterimage" 


6, 7, 8 


3-613.0 12 h 36 m 52.37 s ± 0.06 s 


62°12'25.1" 


± 0.4" 


2.8" 


1.64 


blue irregular 


2, 6, 7 


Other nearby sources within 6": 














3-550.2 12 h 36 m 51.34 s ± 0.06 s 


62°12'26.9" 


± 0.4" 


A R" 
4.0 


1 TO 

1.(2 


blue irregular 


OCT 

2, 0, I 


3-633.1 12 h 36 m 51.60 s ± 0.06 s 


62°12'22.5" 


± 0.4" 


4.2" 


1.72 


irregular ? 


2, 6, 7 


3-659.1 12 h 36 m 51. 72 s ± 0.06 s 


62° 12' 20.2" 


± 0.4" 


5.8" 


0.299 


Sb spiral (Cowie 1999) 


2, 6 


3651+1221 12 h 36 m 51.65 s ± 0.02 s 


62°12'21.4" 


±0.2" 


4.9" 




VLA (20 cm, 3 cm) 


3 


PM3_29 12 h 36 m 51.9 s ± 0.4 s 


62°12'21" 


±3" 


4.7" 




ISOCAM (15 H 


4 


tentative? 12 h 36 m 51.58 s ± 0.06 s 


62° 12' 20.5" 


± 0.4" 


5.9" 




3a contour at 1.3 mm 


1 


Adopted positions for phase calibrators: 






Flux: 


1.3 mm 3.4 mm 




1044+719 10 h 48 m 27.620 s ± 0.002 s 


71°43'35.93" 


±0.01" 


14.0° 




0.39 Jy 0.68 Jy 


5 


1125+596 ll h 28 m 13.342 s ± 0.002 s 


59°25' 14.78" 


±0.01" 


8.8° 




0.12 Jy 0.20 Jy 


5 


1300+580 13 h 02 m 52.465 s ± 0.002 s 


57°48' 37.62" 


±0.01" 


5.5° 




0.13 Jy 0.20 Jy 


5 



Position refs.: (1) This paper; (2) Williams et al. (1996); (3) Richards et al. (1998); (4) Aussel et al. (1999); 

(5) Patnaik et al. (1992). Redshift refs.: (6) Rowan-Robinson (1999); (7) Fernandez-Soto et al. (1999); (8) Zepf et al. (1997) 



tions, then they must be either z < 2 star-forming galaxies 
or, if at higher redshift, heated by active galactic nuclei 
(AGN). To clarify these issues, it is important to obtain 
good positions for the dust sources by millimeter interfer- 
ometry to examine possible coincidences with the Hubble 
Deep Field optical galaxies and VLA radio sources. 

The strongest object in the SCUBA survey is the 
source HDF 850.1, which has a flux density of 7.0+0.4 mJy 
at 850 ^m. Follow-up photometry of HDF 850.1 by Hughes 
et al. (1998) yielded a flux density of 2.1 ± 0.5 mJy at 
a wavelength of 1.35 mm, confirming that the continuum 
is optically thin emission by dust. At longer millime- 
ter wavelengths, the Hubble Deep Field is known to be 
empty of sources. With the BIMA interferometer, Wilncr 
& Wright (1997) found no sources at 2.8 mm, to a 5er- 
limit of 3.5 mJy. Because the 1.3 mm flux density of HDF 
850.1 is close to the minimum detectable level in a 10-hr 
integration at this wavelength, this strongest source in the 
SCUBA survey may be the only source in the Hubble Deep 
Field that can be detected at present in reasonable ob- 
serving times with millimeter interferometers. To measure 
more accurately the position of this source, we observed 
it at 1.3 mm and 3 mm with the IRAM Interferometer on 
Plateau dc Bure, France. 



2. Observations 

The observations were made at 1.3 mm and 3 mm simulta- 
neously, in the interferometer's compact configuration D 
on November 17 and 22, 1998, and in the more extended 
configuration C2 on December 15, 16, and 17, 1998. The 
total integration time was about 40 hours, in excellent 
weather, with precipitable water vapor content s=y 1.5 mm, 
and r.m.s. phase errors at 1.3 mm of < 30°. The five 15 m 
dishes give 10 interferometer baselines from 24 m to 80 m 
in configuration D, and 24 to 183 m in configuration C2. 
At 1.3 mm, the observing frequency was 236.3 GHz, and 
data were taken in upper and lower sidebands separated by 
3 GHz. At 3 mm, data were taken in lower sideband only, 
at 2.8 mm (105.7 GHz) on November 17, and at 3.4 mm 
(88.7 GHz) on all the other dates. The SIS receivers had 
equivalent system temperatures outside the atmosphere 
of 150 K SSB at 3 mm, and 250 to 400 K SSB at 1.3 mm. 
The spectral correlators covered 954kms _1 at 3.4 mm and 
670 km s" 1 at 1.3 mm, with resolutions of 8 and 4kms , 
respectively. Flux densities were calibrated with 3C273 
and MWC349, for which we adopted the following val- 
ues: for 3C273, HJy at 1.3mm and 18 Jy at 3.4mm; for 
MWC349, 1.7 Jy at 1.3 mm and 0.9 Jy at 3.4 mm. The sys- 
tematic uncertainties in the flux scales are estimated to be 
±5% at 3.4 mm and ±10% at 1.3 mm. 
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Fig. 1. IRAM interferometer map at 1.3 mm (236.3 GHz) of 
the continuum emission in a 24" field centered on the SCUBA 
position of the submm source HDF850.1 (cross; the arms of the 
cross correspond to the 3a positional uncertainty of Hughes et 
al. 1998). The synthesized beam is 2.1" x 1.7" (box at lower 
left). The map has not been corrected for the response of the 
primary beam of the 15 m dishes at 1.3 mm (20" FWHM; 
dashed circle). The first positive and negative (dashed) con- 
tours are at ±0.5 mjy (2 a), and the contours thereafter are in 
steps of 0.25 mjy (Iff). 



Fig. 2. Interferometer map at 3.4 mm (88.7 GHz) of the con- 
tinuum emission in a 64" field centered on the original SCUBA 
position (cross) of the submm source HDF850.1. No source is 
detected, to a limit of 0.36 mjy (3a). The synthesized beam is 
4.7" x4.2" (lower left). The map has not been corrected for the 
response of the primary beam at 3.4 mm (53" FWHM; outer 
dashed circle). The inner dashed circle shows the primary beam 
at 1.3mm (see Fig. 1). The first positive and negative (dashed) 
contours are at ±0.24 mjy (2a), and the contours thereafter 
are in steps of 0.12 mjy (la). 



Phases were calibrated with the radio sources 
1044±719, 1125±596, and 1300±580, for which we 
adopted the positions listed in Table 1. All three phase 
calibrators were observed every 20min, at 1.3 mm and 
3 mm simultaneously, and the data were phase calibrated 
with a weighting by the square of their signal-to-noise ra- 
tios and by the inverse square of the system tempera- 
ture. To correct the amplitudes and the phases at both 
3 mm and 1.3 mm for atmospheric seeing, we used the 
1.3 mm continuum total power measurements to compute 
the changes in electrical path length due to short-term 
fluctuations in the atmospheric water vapour content. For 
the 1.3 mm data, we then adopted the phase vs. time cal- 
ibration curve at 3 mm (scaled to 1.3 mm) and on top of 
this, solved for a second-order phase calibration curve to 
fit the phase residuals at 1.3 mm. 

We used both uniform and natural weighting of the u, v 
plane data to make maps. There is not much difference in 
the final result, and the maps shown here (Figs. 1 & 2) 
are the ones made with natural weighting. The 1.3 mm 
map in Fig. 1 includes all the data, in both sidebands 
and in the compact and extended configurations, and has 
been deconvolved with the CLEAN algorithm by Clark 
(1980), with 100 iterations. 



3. Results 

3.1. Detection at 1.3 mm 

A continuum source was well detected at 1.3 mm 
(236.3GHz), in the D configuration on November 17 
and November 22 in the double sideband data, and in 
the upper and lower sidebands separately. All of these 
data subsets agree in the source position and flux den- 
sity, within the errors. On the map that results from 
adding all the data together (Fig. 1), the flux density 
is 2.2 ± 0.3 mjy. This result agrees with the single-dish 
value of 2.1 ± 0.5 mjy measured by Hughes et al. (1998) 
at 1.35 mm. 

The 1.3mm source is at a position 2.4" west and 0.7" 
south of the SCUBA position of the submm source HDF 
850.1 found by Hughes et al. (1998). The SCUBA detec- 
tion was made with the 15 m JCMT single dish, with a 
14.7" beam at 850 /im. Hughes et al. quoted r.m.s. posi- 
tional uncertainties of 0.7" in each coordinate. The inter- 
ferometer position differs from the SCUBA position at the 
3.6 a level, indicating that there may have been an unsus- 
pected systematic error in the SCUBA positions, which in 
turn may affect some of the other submm/optical associ- 
ations suggested by Hughes et al. (1998). 
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Fig. 3. Visibility plot for the continuum source detected at 
1.3 mm. The plot shows the real part of the visibility ampli- 
tude vs. baseline, for u, u-plane data averaged in circular bins 
of 20 m spacing, with la error bars. The data do not have 
a high enough signal-to-noise ratio to yield an angular size. 
The "best-fit" (solid curve) is for a 2.2 mjy gaussian source of 
FWHP diameter 0.7". The dashed curves show the visibility 
functions at the 2 a uncertainty of the fit, for a 1.7" source 
(lower dashed curve) and a point source (upper dashed curve). 



3.2. Size of the 1.3 mm source 

The source is too weak to measure an angular size. In the 
maps from the longer-baseline C2 configuration alone, the 
flux at the source position is 1.3 ± 0.5 mjy, so the source 
may be partially resolved. Our fits to the data in the u, v- 
plane yield a "best fit" FWHP diameter of 0.7" ±0.5". At 
the 2 a level, the visibility data (Fig. 3) agree with sizes 
ranging from 1.7" to a point source. What size should 
we expect? Since we observe mm/submm dust radiation, 
the source cannot be very hot. For example, a circular 
source at z rj 1.7 (see below) u with a dust temperature 
of 35 K, which has a dust opacity of unity at an emitted 
wavelength of 100/mi, would yield a flux of 2.2 mjy at an 
observed wavelength of 1.3 mm if its diameter were 0.3". 
A less opaque source would be even larger. 



3. 3. No detection of 3 mm continuum 

No source was detected at either 2.8 mm (105.7 GHz) or 
at 3.4mm (88.7 GHz), to 3er limits of 0.5 mjy and 0.4 mjy 
respectively (Fig. 2). This is consistent with the 1.3mm 
continuum being optically thin emission by dust. For a 
typical dust continuum spectral index of 3.5, the expected 



We use Ho = 50kms 1 Mpc 1 and qo = 0.5 in this paper. 



fluxes of HDF 850.1 would be only 0.13mJy at 2.8mm and 
0.07 mjy at 3.4 mm. 

3.4. No spectral lines detected 

Although our main goal was the detection of dust emission 
at 1.3 mm, we chose the frequency in the upper sideband, 
237.824 GHz, to coincide with the CO (9-8) line at z = 
3.36, the tentative redshift assigned by Zepf, Moustakas, 
& Davis (1997) to the optical source 3-577.0 located 1.4" 
east of the z f» 1 elliptical galaxy 3-586.0. No line was de- 
tected in the upper sideband, to a limit of 6mJy beam -1 
in channels 32kms _1 wide in the velocity range —300 to 
+367 km s -1 relative to the frequency of the CO (9-8) line 
at this redshift. We also detected no line in the lower side- 
band, 3 GHz lower in frequency. 

At 3 mm, we initially chose 105.7 GHz to search for the 
CO(4-3) line at z = 3.36. At this frequency, no line was 
seen, to a limit of 2mJybeam~ 1 in channels 30kms _1 
wide, in the velocity range ±400 km s -1 around the fre- 
quency of the CO(4-3) line at this redshift. Later on, 
we tuned the 3 mm receivers to 88.738 GHz to search for 
the CO(1-0) line at z = 0.299, the redshift of the bright 
disk galaxy, 3-659.1, located 10" to the southwest of HDF 
850.1. In velocity channels 34kms~ 1 wide, no line was seen 
to a limit of 2mJy beam -1 , anywhere in our 53" primary 
beam, in the velocity range —450 to ±460 km s -1 cen- 
tered on 88.738 GHz. For a galaxy with lines ~ 200kms _1 
wide, this corresponds to a limit of 0.4 Jy kms -1 for 
the CO(1-0) line flux. For comparison, the Milky Way 
has a CO(1-0) line luminosity of 4 x 10 8 Kkms~ 1 pc 2 
inside the solar circle (Rivolo & Solomon 1988), so at 
z « 0.3, the Milky Way would have a CO(1-0) line flux of 
0.05 Jykms" 1 . 

4. Astrometry 

We estimate the r.m.s. positional uncertainty to be 0.3" for 
the centroid of the 1.3 mm dust source. This limit is set by 
the baseline accuracy of the interferometer, the availabil- 
ity of phase calibrators in this part of the sky, the signal to 
noise ratio, and the seeing. Thus far, the best astrometry 
done with the interferometer was for the source W3(OH), 
with phases calibrated on the quasar 0224±671, located 
5° away. In those projects (Wink et al. 1994; Wyrowski et 
al. 1997; 1999), the global statistical positional repeatabil- 
ity was 0.05" r.m.s., including the systematic errors in the 
baseline geometry. To this must be added the 0.01" uncer- 
tainty in the Jodrell Bank-VLA position of the calibrator 
source. For the current Hubble Deep Field project, our 
strongest calibrator is 14° away (Table 1), so the r.m.s. 
error in position would be three times worse, namely 0.2", 
also with Jodrell Bank-VLA calibrator uncertainties of 
0.01" (Patnaik et al. 1992). 

To reduce errors, we chose the three phase calibrators 
at earlier and later right ascensions, and lower and higher 
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Fig. 4. IRAM interferometer map at 1.3 mm (236.3 GHz) of the dust continuum emission superimposed on a greyscale version 
of the BVI image from the Hubble Deep Field. The map is the same as in Fig. 1, but here we only plot the 1.3 mm contours at 
0.75 mjy (3a), 1.25 mjy (5a), and 1.75 mjy (7a). The cross marked SCUBA indicates the position and 3a uncertainty of HDF 
850.1 as given by Hughes et al. (1998). Small crosses with black dots indicate the positions and la uncertainties of the VLA 
sources (Richards 1999) and the large cross marked ISO indicates the position and la uncertainty of the ISOCAM 15/im source 
(Aussel et al. 1999). Optical objects are identified with their photometric redshifts, from the references in Table 1, except for 
the optical source 3-577.0, which has a tentative spectroscopic redshift of 3.36 (Zepf et al. 1997). 
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declinations than the Hubble Deep Field (Table 1). This 
should formally reduce the noise in the baseline errors and 
in our weighted average of calibrators by a factor ~ \/2 
to \/3- To check internal consistency of the reduction soft- 
ware, we made maps of the calibrator sources, and on the 
maps, the position errors were < 0.1". The uncertainties in 
their measured cm-radio positions from the Jodrell-VLA 
survey are ten times lower (Table 1). Two of the calibra- 
tors are weak, ~ O.Uy at 1.3 mm, but their signal-to- noise 
ratio is quite satisfactory at 3 mm. Because we subtract 
the scaled 3 mm phase fluctuations from those at 1.3 mm, 
we can improve the fit of the phase solution at 1.3 mm, but 
there remains a residual instrumental phase error, which 
brings the astrometric part of the position error to 0.2" 
r.m.s. 

In addition to the astrometric error, there is a statis- 
tical error related to the signal-to-noise ratio. For a point 
source, this position error due to noise is 

A9^(B/2)/(S/N) (1) 

where B is the synthesized beam, and S/N is the signal- 
to-noise ratio. For our map in configurations D + C2, the 
beam B is 2.1", and the signal-to-noise ratio is S/N = 7, 
so the position error due to noise is AO = 0.2" r.m.s. The 
convolution of the astrometric error (0.2"), and the noise 
error (0.2") then gives a root sum square error of 0.3". 

5. Source Identification 

The identification of the 1.3 mm dust source with an opti- 
cal galaxy or a cm-radio source in the Hubble Deep Field 
depends strongly on the accuracy of the measured posi- 
tions. The radio-optical registration of the Hubble Deep 
Field is described by Williams ct al. (1996). Their esti- 
mated accuracy of this absolute registration is 0.4". The 
estimated accuracy of the 1.3mm IRAM position is 0.3". 

Radio: The new 1.3 mm position clearly rules out the 
identification with the optical source 3-577.0 proposed by 
Hughes et al. (1998) from the SCUBA data. It also rules 
out the identification, suggested by Richards (1999), of the 
1.3 mm dust source with the VLA source, 3651+1221, lo- 
cated 4.9"±0.4" to the southwest (see further discussion of 
this source in the Appendix). The 1.3mm dust source co- 
incides with the tentative VLA source 3651+1226, which 
is a 4.5ct detection at 3cm (Richards et al. 1998). 

Optical: On the Hubble Deep Field image (Fig. 4), 
the closest object to the 1.3 mm dust source is the arc- like 
feature comprising the optical sources 3-593.0 and 3-593.2. 
From the optical coordinates published by Williams et al. 
(1996), this arc-like feature is 0.8" ± 0.7" southwest of 
the 1.3 mm dust source, where the 0.7" error is the sum 
of the optical and 1.3mm position errors (Table 1). We 
must take the sum of the errors rather than the root sum 
square because the radio-optical registration uncertainty 
is a systematic error rather than a random error. 




2 4 

Redshift 

Fig. 5. Allowed values of the dust temperature and gas 
mass, as functions of redshift, for an unlensed source, that 
are consistent with the observed 1.3 mm flux of 2mJy and the 
450/im/1300/im flux ratio. "Low" redshift (z ~ 1) solutions re- 
quire low temperature (20 K) and high gas mass (6 x 10 10 Mq). 
"High" redshift (z m 3) solutions can have higher dust temper- 
atures (60 K) and lower gas mass (8 x 10 9 Mq). 

The Elliptical: The next nearest optical object is the 
source 3-586.0, located 1.0" ±0.7" east of the 1.3 mm dust 
source. The optical source 3-586.0 is an elliptical galaxy 
at a photometric redshift of 1.0 < z < 1.2 (Hogg et al. 
1996; Lanzetta et al. 1996 (their source number 3-306); 
Mobasher et al. 1996; Gwyn & Hartwick 1996; Sawicki et 
al. 1997; Wang et al. 1998; Fernandez-Soto et al. 1999 - 
their source 303). The optical photometry of this elliptical 
galaxy is tabulated by Fasano et al. (1998), and its spectral 
energy distribution (SED) is given by Cowie (1999). The 
spectrum rises from the optical regime to a maximum near 
2 /im, where it flattens out at a level of ~ 10 /iJy, consistent 
with the non-detection at a level of 23 fiJy at 15 (Oliver 
et al, in prep., reported by Hughes et al. 1998; Aussel et 
al. 1999). 

From the position alone, the dust source could be the 
arc-like feature or the elliptical. We argue, however, that 
the 1.3 mm dust source is probably associated with the 
arc-like feature. Although we cannot associate the dust 
source with this optical object from the astrometry, we 
now present several physical arguments that support this 
association and disfavor the elliptical. 

5.1. Dust flux ratios, dust temperature, and gas mass 

The source HDF 850.1 is not detected at 450 /an to a 
limit of 21 mJy (Hughes et al. 1998), nor is it in the IRAS 
Faint Source Catalog, nor in the ISO images of the Hubble 
Deep Field at 6.7/um and 15 ^m (Aussel et al. 1999). We 
thus have a limit on the 450 /im to 1300 /im flux ratio, 
where the emission is optically thin (the spectral index is 
- 3). This ratio is S , (450/zm)/S'(1300,um) < 10. Let the 
corresponding observed frequencies (667 and 236 GHz) be 
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v 2 and u\, respectively. Because the shape of the Planck 
function is preserved with redshift, the ratio is the same 
at the emitted frequencies, namely 



V2\ (e 



3 / hv 1 (l+z)/kT d 



1) (1 



-t{v 2 {1+z)) 



S{vi) \vi) (e h Ml+z)/kT d - 1) (1 - e -r{ Vl {l+z))^ 



(2) 



where Td is the dust temperature and t(v) is the dust op- 
tical depth. Experience in millimeter detection of dust at 
high rcdshifts during the past decade shows that sources 
are only detected when they contain large amounts of dust 
(<~ 10 8 Mq) — as in the IR ultraluminous galaxies which 
have dust opacities equal to unity near a wavelength of 
100 /im (e.g., Downes, Solomon, & Radford 1993). At this 
wavelength (frequency 3THz), one may approximate the 
dust opacity as t{v) = (^/3)", where v is the emitted fre- 
quency in THz, and n is the dust emissivity index, where 
n « 1.5 to 2. Regardless of whether the dust is optically 
thick or thin, however, at "low" redshifts (z < 1), the ra- 
tio in eq.(2) can be as low as the observed value of < 10 
only with dust temperatures Td < 20 K. (Fig. 5). 

From the mere detection of dust at 1.3 mm with a flux 
density of 2.2 mJy, we can make a useful deduction about 
the gas mass from the formula (cf. Downes et al. 1992) 



M 



1.6 10 6 S„ 



[M ] (l + z) [Jy] 



D, 



[Mpc] 



[K] 



[THz 



-(2+n) 



(3) 



where M is the gas mass, S is the flux density, Dl is 
the luminosity distance, Td is the dust temperature, v cra 
is the emitted frequency, and n is the dust emissivity in- 
dex, which we adopt to be n = 1.5. The main assump- 
tion is that the emission comes from optically thin dust. 
The equation is for a dust mass absorption coefficient of 
0.11 m 2 kg -1 of dust at an emitted wavelength of 1.3 mm, 
as in dense molecular clouds (e.g., Kriigel & Siebenmor- 
gen 1994, their Fig. 12), and a gas-to-dust mass ratio of 
100. The uncertainties are about a factor of three in each 
direction — higher and lower (for an analysis, see Hughes, 
Dunlop, & Rawlings 1997). Higher gas-to-dust ratios mean 
even larger gas masses. 

To have a flux density of 2mJy at 1.3 mm with the 
low dust temperature derived above, a "low" -redshift 
(z ~ 1) source would need a very high molecular gas 
mass, approaching lO n M0. Lower dust temperatures or 
less opaque dust or higher gas-to-dust ratios would all im- 
ply even higher gas masses. Such an enormous quantity of 
dust and molecular gas would have dramatic consequences 
for star formation, and such a large, gas-rich galaxy at low 
redshift would have been easily seen on the HDF images 
or in the follow-up optical spectroscopy of these sources 
(e.g. Zepf et al. 1997). For comparison, a large galaxy like 
NGC 891, at a distance of 10 Mpc, has a total 1.3 mm dust 
continuum flux of 730 mJy, corresponding to a molecular 
gas mass of 1.5 x 10 9 M (Guelin et al. 1993), about the 
same as that of the Milky Way. If NGC 891 or the Milky 



Way were moved to z = 1, at an emission distance (angu- 
lar size distance) of 1.75 Gpc, then the 1.3mm observed 
flux would be only 36 /xJy, that is, 60 times lower than the 
dust flux we detect. 

At higher redshifts, the observed limit on the 450 fim 
to 1300 /xm flux ratio allows the dust temperature to be 
higher. For example, at z = 3, the dust temperature can 
be 60 K, comparable with the blackbody dust temperatures 
deduced for IR ultraluminous galaxies from their IRAS 
fluxes (e.g. Downes & Solomon 1998). The corresponding 
molecular gas masses deduced from the observed 1.3 mm 
flux also become plausible — a few x 10 9 M Q , even in the 
absence of gravitational magnification. Because the line 
of sight to the source passes so close to the elliptical - 
a good gravitational lens — the observed mm and submm 
fluxes may be gravitationally magnified, and the real gas 
mass even lower. 

In summary, our argument, based on the mm and 
submm fluxes alone, that the dust source cannot be a 
"low" redshift (z < 1) object, is as follows: a) the op- 
tical depth is not high; b) the observed flux ratio may be 
obtained at low redshift only if the dust temperature is 
sufficiently low (< 20 K), and c) the observed flux plus a 
low temperature requires an excessive gas mass. 

5.2. Other arguments against identifying HDF 850.1 with 
the elliptical galaxy 3-586.0 

A number of other facts argue against the optical source 
3-586.0, an elliptical galaxy at z m 1, being the optical 
object associated with the 1.3 mm dust source. 

1 ) The gas and dust masses of ellipticals are generally 
too low for 1.3mm dust detections (M gas < 1O 1O M ). 
Even in nearby ellipticals, there are very few detections 
(see, e.g., the review by Knapp 1999). The three galaxies 
with dust detections > 3<r reported by Wiklind & Hcnkcl 
(1995) have gas masses ranging from 1 x 10 5 to 7x 10 8 M . 
The much larger dust and gas mass predicted by the phys- 
ical arguments in the previous section would be totally 
inconsistent with the gas content of an elliptical galaxy 
like 3-586.0. Such galaxies typically have a total dynami- 
cal mass (mostly stars) of 3 x 10 10 M Q within a radius of 
1 kpc, and, with a de Vaucouleurs profile, enclose a total 
mass of 2 x 10 11 M Q at a radius of 7 kpc ( 0.85" at z w 1). 
This argument from the dust and gas mass needed to yield 
2mJy at 1.3 mm is equivalent to the empirical luminosity 
argument by Hughes et al. (1998) — that the mm/sub- 
mm luminosity of the source HDF 850.1 is two orders of 
magnitude too high to be emission by dust in the z m 1 
elliptical galaxy 3-586.0. Therefore the dust source cannot 
be the z = 1 elliptical galaxy. 

2) Dust has not yet been detected at 1.3 mm in any 
"normal " galaxy at z ~ 1. The reasons are firstly, that 
the redshift is not high enough to shift warm dust emission 
into the millimeter band, and secondly, that the emission 
distance (angular size distance) is close to its maximum. 
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A striking exception is the extremely red, z — 1.44 galaxy 
HR 10 (or ERO J164502+4626.4), for which Cimatti et 
al. (1998) and Dey et al. (1999) find mm/submm dust 
fluxes corresponding to an astonishing far-IR luminosity 
of <~ 10 13 L Q , a molecular gas mass of ~ 1O 11 M0, and 
an apparent star formation rate of 1000 to 2000 M Q yr _1 . 
It will be important to investigate whether this object is 
gravitationally lensed or not. 

3) We argued above that if the object is at z w 1, the 
dust must be cold. So far, however, the mm detections of 
dust in high-z objects are mainly of "warm" dust (T > 
30 K), probably circumnuclear regions heated by an AGN 
and/or a nuclear starburst. Therefore it is most unlikely 
that we are detecting the elliptical at z ~ 1. 

6. Discussion 

6.1. The dust source HDF850.1 as an unlensed source 

Hughes et al. (1998) gave a very plausible interpretation 
of HDF850.1 as a high-redshift starburst galaxy, similar 
to the ultraluminous infrared galaxies at z < 0.3. As dis- 
cussed above, from the interferometer position measure- 
ment and the physical arguments about the gas mass, it 
is likely that the dust source HDF850.1 is associated with 
the optical arc- like feature 3-593.0/593.2. The photomet- 
ric rcdshift estimates of the arc are z — 1.73 (Rowan- 
Robinson 1999) and z = 1.76 (Fernandez-Soto et al. 1999). 
In this part of the sky there seems to be a cluster or a 
group of galaxies at photometric redshifts z w 1.7 (the 
objects 3-593.0, 3-633.1, 3-550.2, and 3-613.0). It is con- 
ceivable that the dust source is at this redshift as well, and 
that it is part of the same object as the optical arc- like fea- 
ture. Such a redshift for HDF850.1 would be well within 
the redshift distribution of the identified sources detected 
so far with the SCUBA array in other fields (e.g., Barger 
et al. 1999a,b; Eales et al. 1999; Lilly et al. 1999). 

No optical lines were detected in the arc by Zepf et 
al. (1997). Why not? A simple explanation is that the ob- 
ject is at 1.2 < z < 2.4, consistent with the photometric 
redshifts, so that no strong lines are shifted into the op- 
tical. As shown by the data of Zepf et al. (1997), the arc 
3-593.0 is not likely to be at z > 3, because it is too blue 
— it shows up nicely in the HST F450W band, and it has 
significant flux in the F300W band. It is not a "U-band 
dropout", that is, the hydrogen along the line of sight 
does not greatly decrease the flux in the F300W band, 
as it would if the object were at z > 3, again consistent 
with the photometric redshift of z ~ 1.7. The arc is much 
bluer than, for example, the object 3-577.0, which has a 
photometric redshift of 2.88 (see, e.g, the colors tabulated 
in Zepf et al. 1997). 

In fact, the data on HDF 850.1, if combined with 
the optical fluxes of the optical arc-like feature 3-593.0 
and a redshift of 1.7, agree rather well with the spectral 
energy distribution of the ultraluminous infrared galaxy 




Emitted frequency (Hz) 



Fig. 6. Spectral energy distribution for the dust source HDF 
850.1 and the optical source 3-593.0 if they are at z « 1.7, 
compared with the spectral energy distribution of the ultralu- 
minous IR galaxy VII Zw 31 at z = 0.054. Data for HDF 850.1 
are from this paper and Hughes et al. (1998). For 3-593.0, fluxes 
are from Fernandez-Soto et al. (1999). For the VLA source 
3651+1226, cm-radio fluxes are from Richards et al. (1998) and 
Richards (1999). For VII Zw 31, IR and optical data are from 
Trentham et al. (1999); mm limits are from Downes & Solomon 
(1998), and the 20cm flux is from Condon et al. (1996). The 
solid curve shows a spectrum of the far-IR emitting dust com- 
ponent, from the model described in the text. The dashed line 
is for an assumed synchrotron spectral index a = —0.75 for 
VII Zw 31. 

VII Zwicky 31. This object, which was interpreted by 
Djorgovski et al. (1990) as a merger-induced starburst 
galaxy, was recently selected for study with the HST by 
Trentham, Kormendy, & Sanders (1999) because its "cool" 
spectral energy distribution indicates it is an ultralumi- 
nous galaxy that derives its IR luminosity from an em- 
bedded starburst only, with no "warm" AGN component. 
Figure 6 shows this comparison with VII Zw 31, if HDF 
850.1 is associated with the optical object 3-593.0, that 
has a photometric redshift of z w 1.7. For VII Zw 31, we 
took the IRAS and HST data from Trentham et al. (1999), 
re-scaled for H = SOkms" 1 Mpc^ 1 . The luminosities of 
the far-IR peak of HDF 850.1 and the optical emission of 
3-593.0, calculated for z = 1.7, correspond quite well to 
the observed luminosity of VII Zw 31. The solid curve in 
Fig. 6 shows the predicted spectrum of the cool dust com- 
ponent from the radiative-transfer model for VII Zw 31 
by Downes & Solomon (1998). In that model, the dust 
temperature varies with galactocentric radius. The denser 
dust that is opaque at 100 ^tm has a temperature of 50 K. 
The millimeter dust emission is transparent over the full 
extent of the source, for which the global spatial average 
of the dust temperature is ~ 35 K. 

There is also a good agreement in the size of 3-593.0 
and the dust region in VII Zw 31. On the Hubble Deep 
Field image, the bright part of 3-593.0 has a radius of 
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0.22", or ~ 1.9 kpc at z « 1.7. The total emission along 
the arclike feature, including the faint optical object 3- 
593.2, has an radius of 0.67", or 5.6 kpc. In VII Zw 31, the 
relevant dimensions are nearly the same. The half-power 
radius of the rotating molecular gas disk in VII Zw 31 is 
1.7 kpc, with an outer radius of 5 kpc (Downes & Solomon 
1998), and the ring of star-forming knots in VII Zw 31 has 
a radius of 3.4 kpc (Mazzarella et al. 1999; image in Tren- 
tham et al. 1999). This good match to the star-forming 
region of VII Zw 31 in both size and spectral energy dis- 
tribution, if HDF 850.1 is at z « 1.7 and associated with 
3-593.0, suggests a plausible interpretation of HDF 850.1 
as an ultraluminous IR galaxy with the same intrinsic 
8^m to 1000 /im luminosity as VII Zw 31, namely, Ltr 
= 2 x 1O 12 L . 

Although the best-fit photometric redshift for the opti- 
cal arc-like feature 3-593.0 is with an Scd galaxy template 
at z = 1.73 to 1.76 (Fernandez-Soto et al. 1999; Rowan- 
Robinson 1999), the redshift likelihood function does have 
a much smaller alias as a starburst galaxy at z = 2.6 
(Fernandez-Soto, private communication) to 3.0 (Rowan- 
Robinson 1999). The fit is not as good as for z = 1.73, in 
particular for the J-magnitude, but it is not implausible. 
This is in fact the case discussed by Hughes et al. (1998; 
see their Fig. 5). The spectral energy distribution would 
still match that of VII Zw 31, as in Fig. 6, scaled up by 
a factor of 3.6 (to correct for the greater luminosity dis- 
tance), except in the optical range, where the arc source 
3-593.0 would then be much more blue than VII Zw 31. 
A better fit for both HDF 850.1 and 3-593.0 would be 
obtained with the other two sources measured by Trcn- 
tham et al. (1999), namely, IRAS 12112+0305 and IRAS 
F22491— 1808. Because of the greater luminosity distance 
at z = 3, the 8 /zm to 1000 /*m luminosity would have to 
be L m = 8.2 x 1O 12 L , which would place HDF 850.1 
among the most powerful of the ultraluminous galaxies, 
like Mrk 231 (see, e.g., the sample of Sanders et al. 1988). 

In summary, the arguments in favor of the interpreta- 
tion as an unlensed source at z ~ 1.7 are the fairly good 
match in both size and spectral energy distribution with 
VII Zw 31, the better fit of the optical template with a 
photometric redshift of 1.7 for 3-593.0, the absence of any 
strong lines in the optical spectrum, which is better un- 
derstood if z <~ 1.7, the apparent cluster of objects in the 
range z = 1.6 to 1.8 in this region, the blue color of the 
arc, and the fact that the total luminosity of 2 x 10 12 L 
is not exceptional in the class of ultraluminous galaxies, 
consistent with the absence of any indication of an em- 
bedded optical quasar or Seyfert 1 nucleus (see, e.g., Kim 
& Sanders 1998; Solomon et al. 1997). 

The arguments in favor of an interpretation as an un- 
lensed source at z ~ 3 have been given by Hughes et al. 
(1998) and Carilli & Yun (1999). To these arguments we 
may add the poorer-fitting alias for a photometric red- 
shift of 2.6 or 3.0. The consequence is that the IR lumi- 
nosity would place the source among the most powerful of 



the ultraluminous galaxies, those containing a quasar or a 
Seyfert 1 nucleus. These objects are fairly rare in the exist- 
ing samples of ultraluminous galaxies (see the luminosity 
histogram of Kim & Sanders 1998; their Fig. 3), and usu- 
ally have a small, bright, symmetric nucleus resembling a 
reddened QSO (Surace et al. 1998). There is no evidence 
for such a bright QSO-like nucleus in the optical image 
of 3-593.0, which instead looks like a star-forming region 
extended over several kpc. 

6.2. The dust source HDF850.1 as a lensed source 

Another possibility is that HDF 850.1 is a lensed source. 
The relative abundance of lensed sources in the submm 
is expected to be large, because the submm flux density- 
redshift relations of distant, dusty, star- forming galaxies 
are flat, leading to a larger fraction of brighter images 
in the submm relative to the optical (Blain 1997). The 
fraction of gravitationally lensed sources at the lOmJy 
level at 850 ^m is > 0.01 (Blain 1998, his Fig. 4; Blain et 
al. 1999), orders of magnitude greater than for the corre- 
sponding optical sources. There is therefore a good chance 
that we are detecting a lensed source, especially if the line 
of sight is close to an elliptical galaxy. 

Hogg et al. (1996) originally proposed that the arc-like 
feature 3-593.0 was a source gravitationally lensed by the 
elliptical galaxy at z ~ 1 because of the 1.8" separation 
of the arc from the elliptical 3-586, the mass of the el- 
liptical as deduced from its luminosity, the appearance of 
the arc, concave toward the elliptical, and the fact that 
other gravitational lenses tend to be ellipticals. It is possi- 
ble that the dust source is part of the same object as the 
optical arc-like feature, and lies near a high-magnification 
point of the lens. 

The non-detection of optical lines in the arc is con- 
sistent with 1.2 < z < 2.4 (S. Zepf, private communi- 
cation and Zepf et al. 1997). The low end of this range 
is close to the redshift of the elliptical, giving an unfa- 
vorable geometry for lensing. The high end of the range, 
however, yields nearly the same geometry as for IRAS 
F10214+4724, where the elliptical is at z = 0.9 and the 
source is at z — 2.3. If the optical arc-like feature is at 
its photometric redshift (z ~ 1.7), then the arc source 
and the elliptical are in a favorable configuration for lens- 
ing. The approximate radius of the circle of images may 
be roughly estimated from the standard formula for the 
Einstein radius for an isothermal, spherical lens 



■ Da(LS) 
186kms _ V D A (S) 



(4) 



(e.g., Peacock 1999, his eq. 4.14), where a v is the line of 
sight velocity dispersion of the lens in kms^ 1 , and Da(LS) 
and Da(S) are the angular size distances from the lens to 
the source and from the observer to the source respec- 
tively. 
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Table 2. Comparison of the HDF 850.1/3-586.0/3-593.0 system with other lensed sources 



Source 



— Elliptical galaxy — 
Lens msi4 rnsi4 tuh+k <t v 
redshift r < 0.2" total total 

z lcns (mag) (mag) (mag) (kms -1 ) 



— Lensed feature — 

Source crit. curve msi4 opt./UV 

redshift radii total magnifi- 

z s (arcsec) (mag) cation 



HDF 850.1 system: 

HDF 3-586.0/3-593.0 



1.0 



23.3 



19.5 



< 340 



1.7? 



1.61 



26.1 



Other lenses: 

IRAS F10214+4724 0.9 

HST 14176+5226 0.81 

Cloverleaf ~ 1 ? 



22.8 20.3 18.5 270 2.286 1.2x0.8 20.44 ~ 100 

19.7 260 3.4 2 x 1.5 25.8 

> 24 > 22.5 ~20.5 < 350 2.558 0.7x0.5 17.5 18 to 30 



HDF 3-586.0/3-593.0 data from Cowie (1999); lens model from Hogg et al. (1996). 

IRAS F10214+4724 data and lens model from Eisenhardt et al. (1996); IR data from Evans et al. (1999). 

HST 14176+5226 data from Ratnatunga et al. (1995); Crampton et al. (1996); lens model from Hjorth & Kneib (1999). 

Cloverleaf data and lens model from Kneib et al. (1998a; 1998b), Turnshek et al. (1997), and Chae & Turnshek (1999). 



In the lower-mass, isothermal, elliptical-potential 
model of Hogg et al. (1996), where only 3-593.0 and 
3.577.0 were a candidate lensed pair, the optical arc was 
close to the critical radius of 1.6", from which they de- 
rived a velocity dispersion of > 340 km s -1 if the elliptical 
galaxy is at z = 1. If the dust source is at z ~ 1.7, and 
is a lensed image 0.8" from the elliptical, then the very 
rough estimate from eq.(4) also yields a velocity disper- 
sion of ~ 330 km s -1 . Hogg et al. (1996) tried to explain 
not only the apparent arc 3-593.0, concave toward the el- 
liptical, but also the object 3-577.0 as its counterimage, 
elongated toward the elliptical. The observations of Zcpf 
et al. (1997) made this candidate counterimage unlikely, 
however, because of the tentative detection of a Lya line 
at z = 3.36 in 3-577.0, and the different -B450 — Vsoe colors 
in the arc and the "counterimage" (the arc is bluer than 
its candidate "counterimage" ; for an alternative interpre- 
tation, see Dickinson 1998). 

How could the HDF 850.1 dust source and the optical 
arc-like feature be at different positions? There are sev- 
eral ways this could happen. In the source that is lensed, 
the dust emission may not be in the same place as the 
UV/optical emission. The dust source and the optical arc 
could be different parts of the same galaxy, or a merger of 
two galaxies, with their light paths deflected to different 
spots by the elliptical. They would intersect the ellipti- 
cal's cusp and fold caustics differently, and have different 
magnifications, and different image shapes on the sky. If 
the source is at z ~ 1.7 and lensed, its luminosity would 
still be that of a starburst galaxy. These objects arc often 
highly distorted, with double nuclei, or multiple knots of 
emission — good candidates to yield non-coincident opti- 
cal and mm/submm images if gravitationally lensed. 

Non-coincident images in different wavebands can be 
expected quite naturally, as shown by the comparison with 
other lensed sources in the next section. 



6.3. Comparison with other lensed systems 

More than half of the high-redshift objects in which dust 
and/or CO has been detected at 1.3 mm are known to be 
gravitationally lensed, with magnification factors w 10 for 
the mm dust emission in some of the sources (e.g. IRAS 
F10214+4724, the Cloverleaf quasar, APM 08279+5255, 
MG 0414+0534, and very likely BR 1202-07). Other 
submm-selected sources are also lensed, with lower am- 
plification factors (e.g., SMM J02399-0136 and SMM 
J14011+0252: Ivison et al. 1998; Frayer et al. 1998; 1999). 

Table 2 compares the lower-mass model for the ellipti- 
cal 3-586.0 from Hogg et al. (1996) with similar models for 
IRAS F10214+4724, HST 14176+5226, and the Cloverleaf 
quasar, in which the lensing elliptical galaxies are at com- 
parable redshifts. An important point in these examples 
is that a gravitational lens can image different parts of a 
source to different spots on the sky in different wavebands. 

In IRAS F10214+4724, the combination of an ellip- 
tical at z — 0.9 and a source at z — 2.3 gives strong 
magnification with a 1.18" asymptotic critical radius for 
the UV light rcdshifted to the optical. The optical arc in 
10214+4724 contains two peaks. In the model of Eisen- 
hardt et al. (1996), the brighter east peak is interpreted 
as the blending of two images merging on the critical 
curve, while the west peak is interpreted as the third im- 
age. A possible analogy in the HDF 850.1 system might 
be the brighter optical source 3-593.0 as the two-image 
blend, and the weaker peak 3-593.2 as the third image. 
The 10214+4724 counterimage is 100 times fainter than 
the main arc, and has not been detected in the cm-radio 
or mm-dust emission. The models by Eisenhardt et al. 
(1996) for 10214+4724 predict a 0.4" radial displacement 
in the if-band arc, due to the larger source region at longer 
wavelengths. They note that an 0.5" infrared source would 
be imaged into an elliptical ring with a position angle per- 
pendicular to that of the optical arc, with the center of the 
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ring offset from the elliptical. In fact, the HST NICMOS 
observations of 10214+4724 by Evans et al. (1999) show 
that the peak position of the arc changes even between 
wavelengths of 1.1 and 2.1 /im. 

In HST 14176+5226, the lens redshift, total magni- 
tude of the elliptical galaxy in the HST F814W band, and 
calculated line-of-sight velocity dispersion of the elliptical 
are very similar to those of the elliptical that lenses IRAS 
F10214+4724. In the detailed model by Hjorth & Kncib 
(1999) for the lensing of HST 14176+5226, the magnified 
images in the Einstein cross configuration actually lie out- 
side the critical line, which is an elongated ellipse. 

In the Cloverleaf quasar, the lensing galaxy has re- 
cently been detected in NICMOS images (Kneib et al. 
1998b), and may belong to a galaxy cluster or group at 
z ~ 0.9. The lensing effect may be due to a combination 
of the elliptical and the galaxy group. The different source 
size in different spectral regions leads to quite different im- 
ages. In the Cloverleaf, the four optical images are circular 
spots, while the four millimeter sources are arcs (Kneib et 
al. 1998a). 

6.4- Possible identification of the dust source with the ten- 
tative VLA source 3651+1226 and implications from 
the radio-FIR relation 

From the positions in Table 1 and Fig. 4, the 1.3 mm 

dust source appears to coincide with the weak VLA radio 
source 3651+1226, listed in the Supplementary Catalogue 
of Richards et al. (1998; their Table 5). At a wavelength 
of 3.5 cm, this 4.5<r source has a (revised) sky flux den- 
sity of 7.5 + 2.2 /zJy (E. Richards, private communication). 
From the map by Richards (1999) one may set a 3a up- 
per limit of 23 fi Jy at 20 cm. The dust emission extrap- 
olated from the 2mJy measured at 1.3 mm, would yield 
only 0.02 /iJy at 3.5 cm, so the cm-radio source must be 
mainly synchrotron emission, not dust. Using this VLA 
source in relation with the cm-radio-to-far-IR correlation 
in low-redshift star-forming galaxies, Carilli & Yun (1999) 
concluded that HDF 850.1 could be at a redshift z > 3, 
because the 350/1.4 GHz index af 5 ° is > 1. 

With regard to this submm-to-radio ratio, it is instruc- 
tive to compare the observed flux density of HDF 850.1 
with that expected from the ultraluminous IR galaxy 
Arp 220, if the latter were at high redshift. The result is 
shown in Table 3, which indicates that Arp 220 could not 
be detected much beyond z ~ 1.7 at 20 cm, even with deep 
integrations with the VLA, a fact stressed by Richards 
(1999). An unlensed Arp 220 would not be detectable with 
the IRAM interferometer at 1.3 mm to the detection lim- 
its reported in this paper, in the range 1.7 < z < 4. It 
would also be 3.5 to 6 times weaker at 850 ^m than the 
flux observed by Hughes et al. (1998) for HDF 850.1. Nor 
would it have been detected by the ISOCAM instrument 
at 15 /im, in the same range of redshift. 



Table 3. Expected unlensed flux densities of Arp 220 at high-z 
Arp 220 expected flux at wavelength: 





VLA 


IRAM 


SCUBA 


ISOCAM 


Redshift 


20 cm 


1.3 mm 


850 ^m 


15 /jm 


z 




(mjy) 


(mjy) 


(MJy) 


0.01818 


3 x 10 5 


180 


800 


8 x 10 b 


1.7 


37 


0.3 


1.2 


18 


3.0 


14 


0.5 


2.0 


1.7 


4.0 


8 


0.8 


2.1 


0.8 


HDF 850.1 










observed 


< 23 


2.2 


7.0 


< 23 



Values are scaled to Arp 220's observed spectrum: 

radio: Sopp & Alexander (1991), 

millimeter. Downes & Solomon (1998), 

far-IR: Rigopoulou et al. (1996); Fischer et al. (1998); 

Klaas et al. (1997), near-IR: Soifer et al. (1999). 



As noted by Carilli & Yun (1999), however, the 
submm-to-radio index has a large scatter, so it should be 
regarded more as a redshift indicator rather than a red- 
shift estimator in the sense of the optical photometric red- 
shifts. In particular, for the high-z dust sources detected 
to date, there is not only a large scatter in the 350/1.4 GHz 
index, but also possible contamination of the radio flux by 
emission from an AGN, and modification of the ratio due 
to differential magnification effects of gravitational lens- 
ing. For example, IRAS F10214+4724, at z = 2.3, has an 
index a? 5 4 ° = 0.7, and APM 08279+5255, at higher red- 
shift (z — 3.9) has nearly the same index, a\ 5 2 = 0.8. Both 
of these indices would be compatible with z < 2, given the 
scatter in the index data. 

7. Conclusions 

1) The IRAM interferometer has detected the source HDF 
850.1, the strongest source in the submm survey of the 
Hubble Deep Field by Hughes et al. (1998). We measure 
a flux of 2.2 + 0.3 m Jy, in agreement with the flux density 
measured at 1.35 mm at the James Clerk Maxwell Tele- 
scope. The 1.3 mm flux, together with our non-detections 
at 2.8 and 3.4 mm, and the SCUBA flux at 850 ^m, are all 
consistent with the interpretation by Hughes et al. (1998) 
that the radiation is optically thin emission by dust. 

2) The improved position from the interferometer 
shows the dust source cannot be the optical source 3-577.0, 
as suggested by Hughes et al. (1998). The measured posi- 
tion also rules out the z = 0.299 galaxy 3-659.1 and the 
low surface brightness galaxy 3-633.1 at z = 1.72 sug- 
gested by Richards (1999) as an indcntification for HDF 
850.1. The dust source may coincide with the tentative 
VLA source 3651+1226. 

3) The dust source may be associated with the optical 
arc-like feature 3-593.0/3-593.2, located 1.8" to the south- 
west of the z w 1 elliptical galaxy 3-586.0. The arguments 
in favor of this interpretation are the 0.8" ± 0.7" sepa- 
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ration of the dust source and 3-593 on the sky and the 
temperature-mass argument developed in the text, which 
disfavors an association of the dust source with the ellip- 
tical. If HDF 850.1 is at z w 1.7, the photometric redshift 
of the optical source, then it could be an unlensed galaxy 
with a spectral energy distribution and linear dimensions 
that are similar to those of the z = 0.054 ultraluminous 
starburst galaxy VII Zw 31. Hence, if at z ~ 1.7 and un- 
lensed, HDF 850.1 must be an ultraluminous IR galaxy, 
with an IR luminosity of Lir, — 2 x 10 12 L Q . 

4) Because of the close proximity on the sky of the el- 
liptical 3-586.0, the dust source HDF 850.1 may be gravi- 
tationally lensed. The arguments in favor of this interpre- 
tation are the arclike morphology of the optical source that 
suggests a gravitational mirage; the similarity of the prop- 
erties of the elliptical galaxy to the ones responsible for 
the lensing in IRAS F10214+4724, HST 14176+5226, and 
the Cloverleaf; and the fact that many of the other high- 
redshift objects detected to date at mm/sub-mm wave- 
lengths in dust and CO are magnified by gravitational 
lensing. The arguments against this interpretation are that 
the arc morphology is not conclusive, and that the candi- 
date counterimage, 3-577.0, has a different B — V color, 
and apparently a different redshift. 

The discussion in this paper illustrates the absolute 
necessity of sensitive millimeter/submillimeter interferom- 
eters for identifying sources of dust emission with cosmo- 
logically distant optical galaxies. The next step in this di- 
rection is the construction of a mm/submm interferometer 
with more than an order of magnitude greater sensitivity 
than presently available. 

Acknowledgements. We thank the operators on Plateau de 



There is a 3-sigma contour near the spiral galaxy on 
the IRAM 1.3mm map (Fig. 4), which we list in Ta- 
ble 1 as a questionable source. The 1.3 mm position is 
1" south of the VLA source. The apparent 1.3 mm flux 
density is 0.75 mJy, but this would have to be corrected 
for the primary beam attenuation to ~ 1.2 mJy if the 
source were real. In this case, the emission could not be 
the extrapolation of the synchrotron spectrum from cen- 
timeter wavelengths, because the flux density at 3.4 cm is 
only 16/itJy (Richards et al. 1998). Could it be dust emis- 
sion from the spiral? The problem with dust emission at 
low redshifts (z ~ 0.3) is that most of the energy is at 
shorter wavelengths, not redshifted into the mm /submm 
bands. The dust in the central parts of the Milky Way, 
or in NGC 891, corresponding to a molecular gas mass of 
1.5 x 10 9 M Q , would yield a 1.3mm flux of only 68/zJy 
at z = 0.299, or ~ 20 times lower than the apparent flux 
on the 1.3 mm map. The apparent 3cr signal would thus 
have to come from an extremely gas-rich galaxy, which is 
not in any way indicated by the HST optical image of this 
relatively nearby object. Since there is also a 3cr negative 
peak at a similar radius from the center of the primary 
beam (Fig. 1), the apparent 3ct positive signal may not 
be real. 
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